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ABSTRACT 

We discuss the close link between intermittent events ('quakes') and extremal noise fluctuations which has been 
advocated in recent numerical and theoretical work. From the idea that record-breaking noise fluctuations 
trigger the quakes, an approximate analytical description of non-equilibrium aging as a Poisson process with 
logarithmic time arguments can be derived. Theoretical predictions for measurable statistical properties of 
mesoscopic fluctuations are emphasized, and supporting numerical evidence is included from simulations of 
short-ranged fsing spin-glass models, of the ROM model of vortex dynamics in type II superconductors, and of 
the Tangled Nature model of biological evolution. 
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1. INTRODUCTION 

The term physical aging was coined nearly three decades ago, to describe the slow drift of physical properties of 
polymers and other soft materials. 1 ' 2 The experimental spin-glass investigations of Lundgren and co-workers, 
and Alba and co-workers 6, 7 spurred further interest, and aging phenomena have ever since been investigated by 
a large and growing scientific community. Thanks to these efforts, it is now widely realized that aging is a key 
feature of non-equilibrium dynamics in complex systems, including, e.g. type II superconductors, 8 glasses, 9-12 
granular materials, 13 and soft condensed matter. 14, 15 The presence of common dynamical features across this 
broad field begs the important question of why microscopically very different systems should age in a similar 
way. 

In the following, we briefly review theoretical arguments and experimental and numerical evidence suggesting 
that extremal, record-breaking, fluctuations play a pivotal role in aging after a deep quench. The emphasis is on 
general ideas, and the original papers 16-20 should be consulted for further details. 

2. BACKGROUND 

Complex systems consist of a large number of degrees of freedom, molecules, spins, particles etc., which are 
coupled through a network of mutual interactions and furthermore interact with an external source of white noise, 
e.g., typically, but not exclusively, 17,21,22 a heat bath. The initial quench, e.g. a rapid change of temperature, 
is likely to produce a state of high frustration and strain, a state which is unable to locally fulfill all constraints 
imposed by the mutual interactions. The ensuing aging process tends to relax or optimize the system, resulting 
in a series of metastable, and gradually more stable configurations being visited. Concomitantly, dynamical 
properties as correlation and response function 'stiffen', i.e. their rates of change decrease with the age. 

Being unable to re-equilibrate within realistic time scales, aging materials display a variety of dynamical 
effects, which arc largely independent of microscopic details, but strongly dependent on the system age t w , 
i.e. the time elapsed from the initial quench. Based on the relative values of t w and of the observation time 
t, two dynamical regimes can be identified: For t <C t w , physical averages are approximately constant, and 
two-time correlations are related to the linear response through the fluctuation-dissipation theorem. This 'quasi- 
cquilibrium' regime strongly resembles the time translation invariant dynamics of thermal equilibrium, except 
that many quantities, e.g. relaxation times, may carry a parametric age dependence. For t 3> t w , correlation 




Figure 1. An idealized rendering of the attractors (wedges) visited by an aging system. The arrows represent 'quakes' 
which lead to the new attractors. Among a large number of possibilities (not shown) the marginally more stable attractors 
are typically selected. These quakes entail a large energy change — the vertical displacement of the wedges — and are hence 
de-facto irreversible. By contrast, the degree of stability, which is symbolized by the vertical distance between the bottom 
ad the top of each attractor, only increases by an infinitesimal amount. The layered structure of valleys into valleys likely 
present in complex systems is omitted for typographical simplicity. 



and response functions no longer obey the fluctuation dissipation theorem and physical averages drift at a rate 
which decreases with increasing t w . 

Among the many interesting facets of aging phenomenology, the presence of so-called intermittent events in 
mesoscopic fluctuation spectra 23-26 most directly demonstrates the occurrence of two types of events during the 
unperturbed aging of glassy systems: i) reversible, quasi-stationary fluctuations near a metastable configuration, 
which resemble true equilibrium fluctuations and have a Gaussian probability density (PDF), as in an equilibrium 
state, ii) Jumps, or quakes which lead the system, or in a spatially extended system, parts of it, from one 
attractor, or 'valley' to the next. These sudden and large configurational re-arrangements are dominant in the 
non-equilibrium aging regime, and carry the drift of macroscopic observables. Their fingerprint is a non-Gaussian, 
usually exponential, tail of the PDF of the mesoscopic fluctuations. 

The analysis of several specific examples 16-20 ' 27-29 has shown the intimate connection between quakes and 
extremal, record-breaking, noise fluctuation within metastable attractors. Assuming that extremal fluctuations 
trigger the quakes enables one to translate the properties of record statistics into predictions of the statistical 
properties of the intermittent events, which are found to be in reasonable agreement with experimental and 
numerical results. 18-20 



3. RECORD-BREAKING FLUCTUATIONS AND MEMORY IN GLASSY 

DYNAMICS 

How can record-breaking fluctuations have an impact on the time evolution of physical averages? This would e.g. 
be plainly impossible in thermal equilibrium, where averages do not change. As detailed below, two widespread 
properties of glassy systems, which are sketched in Fig. 1, suffice to produce the effect: i) the presence of numerous 
metastable 'valleys' (defined below), and ii) the 'marginal stability' of these valleys against noise fluctuations. 27 

To fix our terminology, metastable configurations of a noisy dynamical system correspond to stable configura- 
tions of the noiseless dynamics. E.g. for thermal hopping, local energy minima, often called Inherent Structures 
(IS), 30 are metastable configurations. The exit time associated with the basin of attraction of a metastable 
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configuration i.e. the time scale of the escape probability density, provides a quantitative measure of the degree 
of stability, or depth, of the basin. In the case of spatially extended systems with short-ranged interactions, 
similar considerations apply to each of several independently aging subsystems, whose linear size is of the order 
of the thermal correlation length. 

To describe the observed stiffening of dynamical properties in aging systems, it is convenient to distinguish 
between the configurational changes leading to attractors of lesser or equal stability and those leading to attractors 
of greater stability. 16 Changes of the first kind are regarded as internal rearrangements of the 'current' valley, 
while those of the second type, the quakes, lead to a new, deeper, valley. Clearly, then, a configurational change 
qualifying as a quake at the early stages of the dynamical evolution, can at later stages become a reversible 
fluctuation. 

These concepts are translated into a 'non-invasive' exploration algorithm, 31 ' 32 which has uncovered a highly 
non-trivial valley structure of spin-glasses, and shown that the geometrical information culled from time series 
of energy values correlates well with the results of more traditional approaches, using e.g. repeated quenches. 33 
Here, we note that, by construction, valleys grow deeper with increasing t w , and concomitantly include more 
states. They can therefore be expected to acquire a layered substructure of 'valleys within valleys'. 31, 32 This 
hierarchical organization is purely a matter of time scales, as e.g. explicitely shown in an example based on a 
configurational analysis of the Traveling Salesman Problem. 34 

Returning to how record-breaking fluctuations can leave a permanent mark on the evolution of glassy systems, 
we note that condition (i) ensures that a 'sufficiently large' fluctuation is able remove the system from its current 
metastable configuration. The stronger property (ii) means that a fluctuation exceeding all preceding fluctuations 
by an arbitrary small amount, hence a record-breaking fluctuation, is able to elicit a quake. In turn, this implies 
that the rank, not the size, of fluctuations is crucial for aging dynamics. The ability to retain a memory of the 
largest noise fluctuation occurred during the previous evolution can arise through a simple entropic mechanism: 
As the initial quench typically leads to a very shallow valley, we can safely assume that the number of valleys 
of depth s decreases quickly, e.g. exponentially, with s. The marginal increase of the degree of stability of the 
valley selected in a jump must be then understood in a statistical sense, and simply reflects the overwhelming 
predominance of shallow attractors. 16 In this way the size of the largest noise fluctuation experienced so far 
remains imprinted in the height of the barrier of the current attractor, and the temporal sequence of record- 
breaking noise events determines the overall dynamical evolution. 27 

Finally, we note that the importance of marginally stable attractors for complex dynamics was noticed — in a 
more restricted context — by Tang et al. 35 in a seminal paper dealing with a model of elastically coupled degrees 
of freedom, which was devised 36 to describe pulse-duration effects in charge density wave systems. These works 
initiated a decade of investigations in self-organized criticality, 38 ' 39 the study of spatially extended, slowly 
driven complex systems, which spontaneously organize within a set of marginally stable states. SOC focuses on 
the power-laws encountered in the stationary, scale-free, time evolution of systems whose physical properties do 
not change in time, while the log-Poisson theory of aging emphasizes the gradual change of the stability of the 
attractors 29 visited in the aging process. 

4. RECORD-DYNAMICS AS A POISSON PROCESS 

In the simple picture sketched above, quakes are triggered by record-breaking noise fluctuations within each valley. 
A quake entails a large configuration rearrangement, which expectedly dwarfs the extremal fluctuation which 
initiates it. Thus, the most obvious manifestation of record induced dynamics lies not in the size distribution 
of the extremal fluctuations, but in their temporal distribution. Concretely, we seek the time dependence of the 
probability that precisely n noise records — and hence quakes — will occur in the interval (t w , t w +t]. 

To obtain analytical expressions, the correlation between successive noise events must be neglected. This 
approximation is justified when the noise is small and quakes are well separated in time. For thermal noise, our 
description is thus restricted to a range of 'sufficiently low' temperatures. 

With integer times and disregarding correlations, the reversible fluctuations are akin to a sequence of ran- 
dom numbers drawn independently at times 1, 2, 3, . . . , t w , . . . t w + t from the same distribution. Crucially, the 
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Figure 2. The effect of a quake in a two dimensional Ising spin glass model^ is visualized by plotting in white the 
overturned spins. The figure is obtained by projecting the two lowest-lying states of valleys n and n + 1 onto each other. 
The data pertain to a 2d system of size N = 50 2 at T = 0.3, with n = 21. Note the presence of several disconnected 
areas, each comprising a relatively large number of spins. 
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Figure 3. The number of flux vortices entering the system according to the ROM model of flux creep in type II 
superconductors. 19 This model is chosen to illustrate a generic feature of glassy system: the quakes occur in a stepwise 
fashion. Note the logarithmic scale of the abscissa, and the fact that the number of steps per unit of logarithmic time is 
nearly constant.. The age of the system is t w = 100, and the observation time t exceeds a million MC steps. We refer to 
the original article 19 for further details. 
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exact nature of the noise distribution from which the numbers are extracted is immaterial. 22 ' 27 Precisely this 
mathematical property is at the root of the generic features of aging phenomenology. 

The probability that n records occur in a sequence of random numbers between the t w 'th to the t w + t'th 
attempt, is approximately described by the 'log-Poisson' distribution 16, 22, 27 

P n (t w ,t w + t)= ( ttl °g(l+ */*«>))" (1 + t/t w )~ a , t>0; t w >l. (1) 

Strictly speaking, the expression only holds for n small compared to i. 40 We gloss however over this restriction, 
and note that Eq. 1 is a Poisson distribution, whose average 

(n(t w ,t)) = alog(l+ t/t w ) (2) 

depends on the difference of the logarithms of the endpoints of the observation interval (t w ,t w + t], rather than 
on the linear difference t usually encountered. The free parameter a is interpreted as follows: The number of 
records in p = 1,2... sequences updated independently and in parallel is itself a log-Poisson process with a = p, 
and the value a > 1 can cover spatially extended systems, where several regions can evolve in parallel. To include 
the possibility that not every noise record triggers a quake, a is also allowed to assume non-integer values. 

5. MEASURABLE PROPERTIES 

The relevance of record-breaking for aging dynamics has been verified numerically for several models. 17, 19, 27-29 
Below, we mainly discuss predictions which are testable by a statistical analysis of mcsoscopic fluctuations. It 
is assumed that relevant time and age dependent statistical fluctuation properties of a glassy system can be 
obtained through an ensemble of independent measurements, and that these measurement yield a PDF with a 
'normal', typically Gaussian part, which stems from the reversible fluctuations, and an 'anomalous', intermittent 
tail solely due to non-equilibrium intermittent events. Exploiting the large difference in typical size between the 
two types of events, each data point is classified as either a fluctuation or a quake, and the number of 'Gaussian' 
and intermittent events observed are denoted by no and m respectively. 

The shape of the intermittent tail and the relative weight of the two parts clearly depend on the number 
of quakes falling within the observation interval. This influence is exclusively expressed through the average 
(n(t w ,5t)), the only parameter entering Eq. 1. 

We first consider a sufficiently small observation interval St « t w , i.e. such that the St dependence of 
(n(t w ,5t)) can be linearized as 

(n(t w ,St)) = ^. (3) 

Accordingly, PDF data obtained for different observation times St and different ages t w collapse, whenever the 
ratio 5t/t w remains constant. This effect is demonstrated in the lower insert of Fig. 5, which displays the PDF 
of the heat transfer in the three dimensional Edwards-Anderson short-ranged spin-glass model 20 over a time 
interval St for several pairs of values of St and t w with St/t w = 0.01. 

Secondly, we consider varying the observation time t, while keeping t w at a fixed value. The number of 
intermittent events observed grows on average as nj oc log(I + t/t w ). Reversible fluctuation occur at a constant 
rate, and their number grows as no oc t. Hence, the ratio nj/no, which is, basically, the ratio between the 
statistical weight of the tail and the bulk of the PDF, is given by 

nil < X ] 0g (i+ t /t w ). (4) 

n G t w 

The plots in the upper insert of Fig. 5 confirm the linear relationship predicted by Eq.4. The circles are simulation 
data, which are are plotted versus log(l +t/t w ), together with a linear fit (the full line). The main panel in the 
figure shows that the energy e which is liberated in a single quake has an exponential distribution. The above 
properties are all independent of the temperature, as e.g. observed in short-ranged spin glass models. 18, 20, 31 To 
illustrate the independence, we show the rate at which flux vortices enter the system in the Random Occupancy 
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Figure 4. The data plotted describe the distribution of the heat exchanged in spin-glass model 20 over a small (dimen- 
sionless) time interval St, ranging through the values 10, 20, 40, 60 and 80 for the five data set shown. The system has 
age t w = 4000, and the statistics is collected over a time interval t = 2000. The insert show the three central moments, 
average variance and skewness, of the distribution, as a function of St. 




Figure 5. All data plotted describe properties of the PDF of the heat exchange in a spin-glass model at temperature 
T = 0.3. 20 Main panel, and lower insert: The data collapse shows that the statistical properties of the heat exchanged 
in a system of age t w over a time interval St <C t w only depend on the ratio St/t w . In the main panel, the probability of 
energy e being released in a single intermittent event is shown to decay exponentially with e. These data are extracted 
from the tail of the distribution. The upper insert is a numerical verifications of Eq.4. The circles are data points, which 
fall on a fitted straight line, as predicted by the theory. 
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Figure 6. The rate of quakes a(T) as a function of temperature for the ROM model of flux creep in type II supercon- 
ductors. The rate is seen to be nearly temperature independent over a wide range of temperatures. The plot is taken 
from Oliveira et al. 19 



Model (ROM), a thermally activated model of vortex dynamics in type II superconductors 8, 41 which reproduces 
the essential experimental features of these systems. A detailed description is given by Oliveira et al., 19 from 
which Fig. 3 and Fig. 6 are borrowed. The vortex creep occurs in the stepwise fashion shown in Fig. 3, where 
each step corresponds to a quake leading to a new metastable configuration. Fig. 6 shows that the rate at which 
quakes occur remains constant over a wide range of temperatures. 

An activated process with temperature independent rates sounds at first as an oxymoron, but a second look 
gives insight into the subtle role played by the local density of states (LDOS) within each metastable attractor, 
and, more generally, into the memory behavior of complex systems: Importantly, record dynamics only claims 
T independence for the rate at which extremal barriers are climbed. Unlike generic barriers, these barriers are 
determined by the noise history, and hence grow with t w in a temperature dependent fashion. As the temperature 
increases, this barrier growth and the decrease of the Arrhenius time scale for barrier crossing must cancel in 
order to produce the desired effect. Specifically, to get a cancellation, the LDOS for all the attractors must have 
the exponential form 20 

V(e) oc cxp(e/e ). (5) 

Note that thermal metastability is then restricted to the temperature range T < eq. 

The Boltzmann distribution for an exponential LDOS is given by Pe(6) = exp(— ba(T)), where a(T) = 
(1/T — 1/eo) > is an effective inverse temperature. Secondly, with a constant microscopic attempt rate, the 
largest fluctuation on a time scale t w is the largest of 0(t w ) values drawn independently from This extremal 
fluctuation is identical in size to the extremal barrier characterizing the attractors visited at age t w . Applying a 
well known mathematical result, 42 the typical extremal barrier climbed at time t w has height 

b max (T,t w ) oc \n{t w )/a(T). (6) 

The corresponding Arrhenius time is therefore 

t a oc exp(a(T)b max ) = t w , (7) 

as anticipated. 

Interesting effects arise when an aging system is perturbed by small variations of external parameters as 
e.g. the temperature. 18 ' 43 A small negative temperature shift to T' < T imposed at t w does not change the 
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Figure 7. The figure illustrates the effect of a small negative temperature shift (main panel, upper data set) imposed 
on an aging spin glass model 18 and contrasts it to isothermal aging (lower data set). All PDF describe the heat transfer 
over short time St, H(St, t, t w ) = E(t w + (n + l)5t) — E(t w + nSt), n = 0, 1, 2 . . ., where E is the total energy. Main figure, 
lower data sets: The five sets of data shown correspond to ages t w = 100<5t with St = 10, 20, 40 . . . 160 and observation 
window t = t w /2. The temperature is T = 0.20. Upper data sets: One additional cooling step takes the system from 
T = 0.25 to T = 0.20. For the i'th data set, i = 1, 2, ... 5 the negative shift occurs at age t mi , chosen such that the 
corresponding effective ages obtained according to Eq. 8 are t e J f = 2000, 4000, 8000, 16000 and 32000. The corresponding 
sampling intervals are St = t^'/lOO an d * = t%ff/2. Keeping the ratio Stjt e J^ constant collapses the data, whence the 
t e J^ plays the same role as the t w does in the isothermal case. The insert shows that the skewness of the PDF's falls into 
two distinct groups corresponding to T shifted and isothermal aging, with only a small variation within each group. 

extremal barrier already established, but affects the spectrum of thermal fluctuations available to cross it. 18 We 
immediately find 

t a oc exp(a(T')b max ) = t x w , (8) 
where the exponent x = > 1 defines the effective age 

t e J f =t x w , (9) 

The above argument docs not cover a positive temperature shift T' > T. In this case, the ensuing stronger 
fluctuations quickly overcome the extremal barrier b(t w , T) previously established, and reach the extremal barrier 
b(t w ,T') for isothermal relaxation at the final temperature T'. 18 

An approximately exponential LDOS is found in numerical investigations of various complex systems. 33 ' 34, 44 ~ 4f 
Secondly, since the LDOS determines the quasi-equilibrium fluctuation spectrum, it is in principle available from 
calorimetry measurements, which for a spin-glass system lead to the same conclusion. 20 Thirdly, since a quake 
is only irreversible on the time scale at which it occurs, the same process can occur reversibly within an enlarged 
valley at a later stage. Hence, the exponential distribution of the energy released in a quake, see e.g. Fig. 5, is 
itself concurring, albeit indirect, evidence for the ubiquity of exponential density of states in complex systems. 

Let us now consider the statistics of the time intervals between intermittent events. Its precise determination 
requires the ability to discriminate between fluctuations and quakes, and improved experimental techniques are 
required to do so. 47 Within a log-Poisson theoretical description, the theoretical probability to find an interval 
t r between consecutive quakes/intermittent events which is smaller than t is given by 16 ' 31 

R(t\t w ) = l-(l + ^)- a (10) 
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Figure 8. The empirical distribution of the residence time in a valley t r , plotted vs. the scaled variable t/tt^ in the main 
panel and vs. t/t w in the insert. A wide range of t w values was considered. The system in question is a spin glass on a 
fc-regular random graph, but the same type of results is found for Euclidean lattices in 2 and 3 spatial dimensions. 

where a > 1 was introduced in Eq. 1. Note that the average time between the events is (t r ) — t w /{a — 1), which 
is finite for a > 1. Theoretically, the distribution depends on the ratio t/t w only, which is so-called pure aging 
behavior. Equation 10 has been checked numerically 31 with a typical result plotted in Fig. 8: The insert shows 
that a reasonable collapse of the empirical data can be obtained using t/t w scaling. At this level of approximation, 
the analytical expression cannot be distinguished from the data. The better collapse shown in the main panel 
is obtained with the scaling variable t/t 1 ^, where fi is a weakly size dependent parameter close to unity. For the 
chosen parameters, /i = 1.055. Other numerical simulations 48 indicate that the small deviations from t/t w scaling 
decrease systematically as the system size grows and/or the temperature decreases. Experimental deviations are 
observed in several other cases. 15 In spin-glasses they are attributed to the finiteness of the cooling rate. 49,50 
Such effects are not covered by this theory in its present form. Note that Figs. 2 and 8 are not obtained by an 
intermittency analysis, but stem from the numerical exploration algorithm 31 previously mentioned. The close 
relation between the two approaches is briefly discussed below. 

As further evidence of the relevance of the Log-Poisson description for complex dynamics, we finally consider 
the distribution of the 'log waiting times', which are defined as follows: Let t\ . . . tk be the times at which quakes 
1 . . . k occur. The waiting times between the events are independent and exponentially distributed stochastic 
variables in a usual Poisson process, while in a Log-Poisson process the same is true for the 'logarithmic waiting 
times' defined as 

r k = log(i fe ) - log(i fe _i) = log(i fe /^_i). (11) 

In principle, this distribution is measurable in intermittency experiments by recording the intervals between 
anomalous spikes. Numerical results for the ROM model previously described and for the Tangled Nature model 
of biological evolution 51 are shown in Fig. 9. 

6. SUMMARY AND CONCLUSION 

The Log-Poisson non-equilibrium aging scenario leads, in a general and near parameter-free manner, to analytical 
expressions for several statistical properties of mesoscopic fluctuations which are of relevance for the intermit- 
tent behavior of glassy systems after a deep quench. The distinction between time homogeneous fluctuations 
and time heterogeneous quakes is built right into the description, which hence naturally covers various types of 
complex dynamics, where the same distinction applies, e.g. driven dissipative systems, 27, 29 and biological evo- 
lution. 17 ' 21 ' 28 ' 52 ' 53 Broadly speaking, the physical mechanism behind all these examples combines an entropic 
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Figure 9. The empirical distribution for the log waiting times numerically obtained for the ROM and the Tangled Nature 
model 17 ' 51 is shown. In good agreement with Eq.ll, the results nearly follow an exponential decay (broken line). 

aspect, the over-abundance of attractors with low stability, and a dissipation (of e.g. energy) mechanism, which 
guarantees the ability of the system to spontaneously drift away from its initial state. These aspects have not 
received much attention in other theories of aging, where the main focus is on the linear response to a small 
external perturbation. 54 ' 55 In spite of this important difference, the log-Poisson description also adapts and/or 
incorporates other established ideas about complex dynamics. The emphasis on marginal stability modifies 
and extends pre-SOC developments 35, 36 on memory behavior. The increasing level of hierarchical organization 
acquired by the valleys is a recurring theme in complex dynamics. 56, 57 Finally, attractors (and valleys) are 
defined at the level of thermally correlated clusters. 20 These are localized objects separated by frozen degrees of 
freedom, and their increasing stability arguably corresponds to the growth of the thermal correlation length. 12 

Several questions are under investigation, e.g. more detailed predictions for the fluctuation spectrum of 
correlation estimators used in measurements, 15 the origin of the exponential density of states and extending the 
picture to the case where attractor properties beyond the stability can change as the system ages. 
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